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Abstract 


We present an approximation scheme to solve the Non Perturbative Renormalization Group 
equations and obtain the full momentum dependence of the n-point functions. It is based on an 
iterative procedure where, in a first step, an initial ansatz for the n-point functions is constructed 
by solving approximate flow equations derived from well motivated approximations. These ap¬ 
proximations exploit the derivative expansion and the decoupling of high momentum modes. The 
method is applied to the O(N) model. In leading order, the self energy is already accurate both 
in the perturbative and the scaling regimes. A stringent test is provided by the calculation of 
the shift A T c in the transition temperature of the weakly repulsive Bose gas, a quantity which is 
particularly sensitive to all momentum scales. The leading order result is in agreement with lattice 
calculations, albeit with a theoretical uncertainty of about 25%. 
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I. INTRODUCTION 


The need for reliable and efficient non-perturbative calculation methods is felt, in various 
forms, in nearly all fields of physics: in nuclear and particle physics (to deal with the infrared 
sector of quantum chromodynamics and the associated phenomena of color confinement and 
chiral symmetry breaking), in condensed matter and statistical physics of systems in or out of 
equilibrium (phase transitions and critical phenomena, disorder systems, strongly correlated 
electrons), to quote just but a few general examples. In many of these cases, the absence 
of a small parameter prevents one to build a solution in terms of a systematic expansion. 
In order to treat such problems, what one needs is a non-perturbative method that allows 
the calculation of correlation functions for arbitrary values of the external momenta, from 
which most physical quantities can be deduced. 

Among the non perturbative methods that have been developed along the years, the non 


.1, y, k y, q 


stands out as a very promising 


perturbative renormalization group (NPRG) 
tool, suggesting new approximation schemes which are not easily formulated in other, more 
conventional, approaches in field theory or many body physics. The NPRG has been ap¬ 


plied successfully to a variety of p 
physics (for reviews, see e.g. j(J 
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'hysical problems, in condensed matter, particle or nuclear 
). In most of these problems however, the focus is on 
long wavelength modes and the solution of the NPRG equations involves generally a deriva¬ 
tive expansion which only allows for the determination of the n-point functions and their 
derivatives at small external momenta (vanishing momenta in the case of critical phenom¬ 
ena). In many situations, this is not enough: a full knowledge of the momentum dependence 
of the correlation functions is needed to calculate the quantities of physical interest (e.g. to 
get the spectrum of excitations, the shape of the Fermi surface, the scattering matrix, etc.). 

The NPRG presents itself as an infinite hierarchy of equations relating sequentially the 
various n-point functions. To our knowledge, most efforts to solve this hicrachy, aside 
from the derivative expansion alluded to above, have been based on various forms of the 


early proposal by Weinberg Q, that is they involve some truncation of the infinite tower 
of flow equations for the n-point functions, ignoring higher order vertices, or possibly using 
various ansatzs for some of them fill fill ]. This leads to approximations similar to those 


used when solving the hierarchy of Schwinger-Dyson equations |12j. However, despite the 
fact that very encouraging results have been obtained this approximation scheme presents 
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convergence dificnlties 113]. 

The goal of this paper is then to present a method for solving the NPRG equations that 
keeps the contribution of all the vertices present in the flow equations. This is achieved by 


exploiting specific properties of the NPRG. The method allows one to get, in a relatively 


simple way, the full momentum dependence of the n-point functions. It involves iterations 
that starts with an initial guess for the n-point functions. That initial guess is then injected 
in the flow equations which are integrated in order to obtain a leading order expression 
for the n-point functions. And so on. Clearly, each new iteration involves more n-point 
functions, and the scheme may become rapidly prohibitively complicated. It is therefore 
crucial that the starting point of the iterations, that is, the initial ansatz for the n-point 
functions, be as close as possible to the exact solution, in order to get a good approximation 
with a minimum number of iterations. The construction of this initial ansatz is therefore 
the central part of the method. 

To derive this initial ansatz we shall first simplify the flow equations using well motivated 
approximations. We shall exploit a modified derivative expansion in its leading order and 
the decoupling of high momentum modes in the flow equations in order to simplify the 
momentum dependence of the vertices that govern the flow. The resulting approximate 
equations are then solved explicitly. 

The particular class of problems that we are concerned with can be formulated in terms of 
a field theory, and as a generic case, we shall consider here a scalar </> 4 theory in d dimension 
with O(N) symmetry: 



( 1 ) 


where the field ip(x) has N real components with i — 1, ■ • • , N. We emphasize however 

that most of the arguments presented in this paper have a wider range of applicability. 

In this paper, we shall apply the method to the calculation of the self-energy at criticality 
and at zero external field, in leading order and in d — 3. This involves getting the initial 
ansatz for both the propagator and the 4-point function. Constructing this initial ansatz 
is the main task carried out in the present paper. It is presented in sect. ITTTI together 
with a more detailed description of the approximation scheme. First, in sect. m we review 
basic features of the NPRG, and illustrate various strategies that have been used to obtain 
solutions of the flow equations. These will provide the necessary background to motivate 
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the approximation scheme presented in sect. nn as well as the approximations involved in 
the construction of the initial ansatz for the 4-point function. The reader familiar with the 
NPRG may skip this section and go directly to sect. 11111 The results for the self-energy are 
presented in sect. I1VI 

The self-energy thus obtained has the correct behavior at all momenta. It agrees with 
perturbation theory in the ultraviolet and it presents the expected power-law behavior in 
the infrared. As a benchmark for our approximation scheme we shall use the shift A T c 


(see also 


for a 


of the transition temperature of a weakly interacting Bose gas 00 
recent review on the theory of the weakly interacting Bose gas). As we shall recall later, 
the precise evaluation of A T c requires an accurate knowledge of a 2-point function at all 
momentum scales, and it constitutes therefore a very stringent test of any methoc 
at getting the full momentum dependence of n-point functions. As shown in Ref. 
calculation of A T c reduces to that of the change A(</? 2 ) of the magnitude of the fluctuations 
of the field described by the action o, for d = 3 and N = 2 |12|. This calculation can be 
done immediately once the self-energy is known. It is presented in sect. IIV Bl together with 
a comparison with estimates of this quantity using different techniques, for instance lattice 

3 , 0 . 


aiming 
14], the 


calculations 


In a companion paper 0], we extend the method described here to the next-to-leading 


order calculation of the self-energy (which involves the leading order calculation of the 


However, since the publication of ref. 0], we have been able to improve the accuracy of the 
leading order calculation of the 4-point function, which yields a considerable improvement of 
the next-to-leading order self-energy; the final results that we obtain for A T c are in excellent 
agreement with the lattice calculations, with a much reduced theoretical uncertainty as 


point function). Some of the results of this study have already been presented in ref. 0 


compared with the estimates presented in the present paper Further progress has been 
achieved in an effort to get rid of some of the approximations used in the present work, and 
which contributes to the theoretical uncertainty in the predictions. A possible strategy to 
do so has been presented in 22], and first results concerning its numerical implementation 
will be presented shortly Q- 
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II. SOME FEATURES OF THE NPRG EQUATIONS 


Generalities 


The NPRG allows the construction of a set of effective actions T K [0] which interpolate 
between the classical action S and the full effective action P[d]: In r,.[</>] the magnitude of 
long wavelength fluctuations of the field is controlled by an infrared regulator depending on 
a continuous parameter k which has the dimension of a momentum. The full effective action 
is obtained for the value k — 0, the situation with no infrared cut-off and where therefore 
all fluctuations are taken into account. In the other limit, corresponding to a value of k of 
the order of a microscopic scale A at which fluctuations are suppressed, T K=A [0] reduces to 


the classical action 


421- 


In practice the control of the magnitude of the fluctuations is implemented by adding to 
the classical action m the regulator 


j 


d d q 
(2tt) c 


ipi(q) R K (q) tpi(-q), 


( 2 ) 


where R K denotes a family of “cut-off functions” depending on k. The role of A S K is to 
suppress the fluctuations with momenta q ^ k, while leaving unaffected those with q > n. 
Thus, typically, R K (q) —> n 2 when q <C n, and R K (q) —■> 0 when q ^ n. There is a large 


freedom in the choice of R K (q), abundantly discussed in the literature (24 l25l . 126. £]. The 
choice of the cut-off function matters when approximations are done, as is the case in all 
situations of practical interest. We have used in this work the cut-off function proposed in 


24]: 


RM 2 )^{K 2 -q 2 )e{K 2 -q 2 ). 


(3) 


This regulator allows many calculations to be done analytically. It is known to work well 
with the derivative expansion in leading order, which we shall use in this work. 

For each value of k, one defines the generating functional of connected Green’s functions 


W K [J\ = In / Dip exp < —S^] - A5V,[y>] + / d d x(p(x)J(x) > . 


We have, for instance, 


4>kA x ) = MuiLj = 


dW K 

5J(x) 


(4) 


(5) 
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FIG. 1: Diagrammatic illustration of the r.h.s. of the flow equation of the effective action, eq. ©• 
The crossed circle represents an insertion of d K R K , and the thick line a full propagator in an 
arbitrary background field. 


The Feynman diagrams contributing to W K are those of ordinary perturbation theory, except 
that the propagators contain the infrared regulator. We also define the effective action, 
through a modified Legendre transform that includes the explicit subtraction of A S K : 


U[^] = + / d d x <j>(x)J^,(x) - AS k [0], 


( 6 ) 


where is obtained by inverting eq. ©• Note that, in this inversion, 0 is considered as a 
given variable, so that becomes implicitly dependent on k. 

One can write an exact flow equation for T K [0] which gives its variation as a function of 
k,, at fixed 0. It reads 


3* r «[0] = 7T tr [ TTWi d * R M 2 ) [ r L 2) + R *] a - a > 


(2tt) 


(7) 


where rl 2 ' is the second derivative of T K with respect to 0, and the trace tr runs over the 0(A) 
indices. Eq. © is the master equation of the NPRG. Its right hand side has the structure 
of a one loop integral, with one insertion of d K R K (q 2 ) (see fig. ©. The solution of eq. © 
yields the effective action T[0] = T K= o[0] starting with the initial condition T k= a[0] = S[(p] 
(see e.g. jfj). 

As well known (see e.g. j3|), the effective action T[0] is the generating functional of 
the one-particle irreducible n-point functions. This property extends trivially to T K [0]. 
Since we shall be concerned only with n-point functions for constant (eventually vanish¬ 
ing) external fields we exploit translational invariance to define reduced ?r-point functions 
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FIG. 2: Diagrammatic illustration of the r.h.s. of the flow equation for the 2-point function, 
eq. ©. The black dot denotes the four-point function and the thick line the full propagator G. 
The circled cross represents the insertion of d K R K . 


r^(«;pi,... ,p n ) as follows : 

(2n) d S^ d \p 1 + ---+p n )r^\ K ;p 1 ,...,p n ) 

= [ d d x x ... [ d d .r„ Pj ' Xj 


8 n T h 


S(j>(xi). ..8(j){x n ) 


( 8 ) 


t>=cst 


By deriving eq. (ED with respect to 0, and then letting the field be zero, one gets the flow 
equations for all n-point functions in a vanishing background held </>. For example, the 
equation for the 2-point function reads: 


<9 K r!?( k;p) = 8 l2 djl( k\p) = ~\j jyyy 


d K R K (q) G' 2 (k-, q) r^ z (/c;p, -p, q, -q), 


where we have introduced the self-energy £(k; q) and 

G _1 (/y q) = q 2 + R K (q) + £(k; g). 


(9) 


( 10 ) 


In eq. ©, and later in this paper, we often denote simply by numbers 1,2, etc., the O(N) 
indices ii,i 2 , etc., in order to alleviate the notation. A diagrammatic illustration of the 
right hand side of eq. © is given in fig. [21 Similarly, the how of the the 4-point function in 
vanishing held reads: 


d K T[% i (K;p 1 ,p 2 ,p 3 ,P 4 ) = I yd K R k {q 2 )G 2 (K;q) 

x |g(k; q + pi + p 2 )r^( k; Pi,P2, Q, ~q-Pi~ P2)^u ij (k; Ps, Pa, ~q, q ~ Ps ~ Pa) 
+ G(k; q + Pi+ p 3 )rS^ («; Pi, p 3 , q, -q ~ Pi ~ p 2 , p 4 , -q, q ~ P2 ~ Pa) 


+ G(k; q + pi+ PA^uij («; Pi, p 4 , g, ~q~ Pi- PaF;^^: p 3 , p 2 , -g, g - p 3 - 


(4) 


1 

2 


d d q 5 

■^■p< 9 KJ R K (g)G 2 (/q g)r 4 2 34 j i («;p 1 ,P2,p 3 ,p 4 , g, -g). 


( 11 ) 
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FIG. 3: Diagrammatic illustration of the r.h.s. of the flow equation for the 4-point function, 
eq. EH): contribution of the 4-point functions (represented by black disks) in the three channels 
s, t and u, from left to right. The crossed circle represents an insertion of d K R K , and the thick line 
a full propagator. 

The four contributions in the r.h.s. of eq. m are represented in the diagrams shown in 
figs. 01 and |U 

Eqs. © and (HU> for the 2- and 4-point functions constitute the beginning of an infinite 
hierarchy of exact equations for the n-point functions, with the flow equation for the n-point 
function involving all the m-point functions up to m = n + 2. Clearly, solving this hierarchy 
requires approximations. In the rest of this section we discuss various approximations that 
are commonly used in the context of the NPRG, and that we shall exploit in the more general 
scheme presented in the next section. In the next sub-section we recall how perturbation 
theory can be recovered from the hierarchy through an iterative procedure. Then, we focus 
on the regime of small momenta where an expansion in powers of gradients of the field often 
yield accurate results. In particular we briefly discuss the leading order of this expansion, 
the Local Potential Approximation (LPA). Finally, in the last subsection, we review simple 
properties of correlations functions of the O(N) model in the limit of large N: this will 
provide a simple, yet non trivial, example in which the momentum dependence of correlation 
functions can be analyzed in detail. 
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FIG. 4: Diagrammatic illustration of the r.h.s. of the flow equation for the 4-point function, 
eq. m- contribution of the 6-point function (represented by a black disk). The crossed circle 
represents an insertion of d K R K , and the thick line a full propagator. 

B. Perturbation Theory 


Perturbation theory can be recovered by solving the exact flow equations iteratively, 
starting with the classical action as initial input (see e.g. IsfJ). The perturbative 
expansion of the effective action, or equivalently the loop expansion, is controlled by a 
“small parameter ”, namely h. Making this parameter explicit one rewrites eq. © as 

d K r K [<P] = 5 Tr d K R K [ri 2) + r k ] _1 , (12) 

and then proceed to the expansion in powers of h. In leading order P re [<fi] is independent of 
k and is equal to the classical action: 

rL 0| M = sW + offi). (13) 


The solution of this equation can be used to obtain an approximation for the 2-point function 
T^ [</>], by taking the second derivative of S[(j)] with respect to </> (see eq. ©). At the next 
iteration this is used in the right hand side of eq. m in order to obtain the order one 
correction to T K [</>]. One gets then, after integrating the flow equation from A to n (using 
the fact that S^ is independent of k) 


TW[0] = S[0] + ^ Trlog 


S& + R k 


+ 0(h 2 ), 


(14) 


S^ + R a _ 

where one recognizes the familiar one-loop expression of the effective action. One can repeat 
the procedure and show that, after n iterations, one reproduces the result that one would 
obtain by calculating T K [</>] using perturbation theory at order n-loop (with the IR cut-off). 
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In the case of massless theories, which we are interested in here, this iteration scheme 
is applicable only for values of k not too small. Indeed, in general, perturbation theory 


stops to make sense 


2.8] when k ^ k c , with k c ~ w 1 ^ 4 ^ , where u is the coupling constant 


defined in eq. (JTJ) (concrete estimates of k c will be presented in the next subsection). When 
k —» 0, perturbative calculations may lead to infrared divergent expressions. This difficulty 
is particularly important in the scaling regime, where p -C p c ~ k c . Other approximation 
schemes are then required. 


C. The local potential approximation 

The derivative expansion offers the possibility to calculate some properties of the scaling 
regime. It exploits the fact that the shape of the regulator in the flow equations (e.g. 
eqs. © or m forces the loop momentum q to be smaller than k, i.e., only momenta q < k 
contribute to the flow. Besides, in general, the regulator insures that, as long as hi ^ 0, 


all vertices are smooth functions of momenta 


43]. Then, in the calculation of the n-point 


functions at small external momenta pi, it is possible to expand the n-point functions in 
the r.li.s. of the flow equations in terms of q 2 /hi 2 and p 2 / n 2 , or equivalently in terms of 
the derivatives of the held. Note, however, that since eventually k —> 0, such an expansion 
strictly makes sense only for p t = 0, unless there is a mass in the problem. 

In leading order, this procedure reduces to the so-called local potential approximation 
(LPA), which assumes that the effective action has the form: 


^LPA r 


[</>] = J d d x + V K {p)\, (15) 

where p = 0j0j/2. The derivative term here is simply the one appearing in the classical 
action, and V K (p) is the effective potential. The exact how equation for V K is easily obtained 
by assuming that the held 0 is constant in eq. ©• One needs however to take into account 
the O(N) symmetry, and to decompose the propagator of the scalar held in a constant 
background 0« into its transverse (Gt) and longitudinal (Gl) components: 


Gij(K-,q) = G t (k; q) ^ 

Then the equation for the potential reads: 


2 P 


+ Gl(k ; q)- 


2 P 


d K V K (p) = 


d d q 
(27r) c 


d K R K (q ) {(N - 1 )G t (k] q) + G l (k ; q)} . 


(16) 


(17) 
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By using the LPA effective action, eq. (uni, one gets 

1 


Gt(k',q) — 
G L {iv,q) = 


q 2 + V'(p) + R k (qy 
1 


q 2 + V'(p) + 2pV"(p)+R K (q )’ (1§) 

with V'(p) = dV/dp and V"(p) = d 2 V/dp 2 . With these propagators, eq. (1T7I) becomes then 
a closed equation. 

Higher order corrections to the LPA include terms in the effective action with an increas¬ 
ing number of derivatives. Although there is no formal proof of convergence, the derivative 
expansion exhibits quick apparent convergence if the regulator R K (q) is appropriately cho¬ 


sen |241271 131|. In practice, the LPA reproduces well the physical quantities dominated by 
small momenta (such as the effective potential or critical exponents) in all theories where 


improvements j^], and the derivative expansion has been pushed up to third order 
yielding critical exponents in the Ising universality class of the same level of precision as 
those obtained with the best accepted methods (see e.g. 

An interesting improvement of the LPA, which we refer to as the LPA’, takes into account 
a running field renormalisation constant Z K and allows for a non trivial anomalous dimension, 


it has been tested (see, for example, p, 13 ]). Higher order corrections lead to significant 


Ei, 


determined from the cut-off dependence of Z K 
to be of the form: 

t^LPA! r •’ f •-> ' Z 


Q. 


In the LPA’, the effective action is assumed 


[c>] = / cfx\ + V K (p) 


(19) 


where Z K is a function of k (and not of p). It is useful to explicitly include the held 
normalisation in the regulator ©, he., we redefine R K by multiplying it by the factor Z K . 
Thus, the regulator used in the present work becomes: 


RM = Z« (k 2 - 9 2 )0(k 2 - 9 2 )- 


( 20 ) 


The factor Z K is determined from the how equation for in a constant external held, 
which can be derived from eq. ©• The vertices and propagators entering this equation are 
those dictated by the form (USD assumed for the effective action. By expanding the resulting 
equation to order p 2 \ r^(/c;p) — r^(/q 0) ~ p 2 Z K (recall that for non vanishing k, T^ 2 \n',p) 
is a smooth function of p), one obtains the following equation for Z K 


d n Z K = -p 0 (W / (p 0 )) 2 4 


d d q 

(27r) fl 


q 2 Gl^q)G 2 T ^q)(Z K + R' K (q)y, 


( 21 ) 
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where R' K (q) = dR K (q)/dq 2 , the derivative d K acts only on the explicit factors R K (and their 
derivatives), and p is fixed at its running minimum p = p 0 (which depends on k). The 
anomalous dimension is related to Z K by (see e.g. j(|; for a simple proof, see app. f~Yl) : 

rj K = -Kd K ln Z K . (22) 

In the LPA’ the flow equation for the effective potential is the same as in the LPA, eq. <HZD, 
except for the replacement q 2 —> Z K q 2 in the propagators. It follows that the flow equation 
for the potential is coupled with the flow equation for Z K , eq. (EH). 

The derivatives of V K (p) with respect to p give the n-point functions at zero external 
momenta as a function of k. We shall be mostly concerned in this paper with the critical 
regime where p 0 ( K = 0) = 0, and hence in n-point functions in vanishing external field, for 
which we shall introduce special notation. We set: 




d 2 V K 


; d 3 V K 


m K = 

K dp 

1 

p= o 

Ill 

a- 

to . 

1 

P= 0 

dp 3 

P= o 


For vanishing external field the propagator is diagonal, G i 2 (/c; q) = ^^Glpa 1 ^] q), with 

G LPA'( K ’> s) = Z ^ 2 + r k(q) + m l • (24) 

For the n-point functions and we have, respectively: 

ri.234 (k) = 9 k (A 12 ^34 + <5l3 <^24 + <5l4 $23) j (25) 


and 

41 123456 ( K ) = [^56 (<5l2^34 + <5l3^24 + ^14^23) + <^46 (<^12^35 + <^13^25 + <^23^15) 

+^36 (<^12^45 + ^14^25 + ^15^24) + <^26 (^13^45 + ^14^35 + ^15^34) 

T5l6 (^23^45 + <524(535 + ^25(534)] • (26) 

In order to factor out the large variations of the effective potential which arise when k 
varies from the microscopic scale A to the physical scale k — 0, and also to exhibit the fixed 
point structure, it is convenient to isolate the explicit scale factors (V K ~ K d , Z K p ~ fiW 2 ) 
and to define dimensionless quantities: 

v K (z) = K^K~ d V K (p), (27) 
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with 


z = />;, 1 /, K 2 - d P . 


(28) 


In these definitions, for further simplifications, we have also included a factor K c i , which 
originates from angular integrations: 


K- 1 = 2 d ~ 1 TT d / 2 d r(d/ 2 ). 


(29) 


Note that can be a small number, e.g. AT 3 = 1/67T 2 . We also introduce dimensionless 
couplings: 


ml = Z k k 2 ml g K = K^Z^ a g K , h K = K?Z*k»- m h K , 


— TS — lr72,A—d 


_ —2 v3 ,3—2d 


(30) 


so that: 


rn't: = 


dv K 

d 2 v K 


f d 3 v K 

dz 

9k ~ dz 2 

z= 0 

z= 0 

K ~ dz 3 


(31) 


z=0 


The solution of the LPA’ is well documented in the literature (see e.g. (6, |27|). It is 
convenient to solve the equation for the derivative of the potential with respect to z, i.e., 
w K (z) = d z v K (z), rather than that for the effective potential itself. With the Litirn regulator 
PI). the integrals in eqs. m and m can be done analytically. One gets: 


Kd K w K = -(2-r] K )w K + (d-2+rj K )zw' K - ( 1 


and 




(N—1)w' k 3w' k + 2ziu 


d + 2j\(l + w K ) 2 {l-\-w K -\-2zw' K ) 2 j'' 


(32) 


= 


4-o«(-o)) 2 


(33) 


(1 + 2z 0 w' k (z 0 )) 2 ’ 

where w' K = d z w K (z), w" = d 2 w K (z ), and z 0 = z 0 (fv) is the running minimum of the potential 
(w k (zq) = 0). Eqs. (EH and (ESj) are solved starting from the initial condition at k = A: 


w K (z, k = A) = rh\ + g A z, 

where m A and g A are related to the parameters r and u of the classical action 0 by 

. o r , u K d 


(34) 


m A = ^ 


9a = 


A 4_d 3 


(35) 
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Before looking at some results obtained by solving numerically eqs. (TTZl) and ED, it 
is useful to get insight into the expected behaviour of the solution by solving eq. (H?!) 
approximately JfJ, ignoring the anomalous dimension. To this aim, we assume that, for all 
k, w K (z) retains the form of eq. (ED, he., 

w K (z) = ml + g K z. (36) 


The minimum of the potential, z 0 (k;), satisfies w K (z 0 ) = 0, i.e., z 0 (k) = —m 2 K /g K . The 
equations for z 0 (k ) and g K are easily obtained from eq. ED, taking into account that, at 
criticality, z$g -C 1 to make simplifications whenever appropriate. One gets: 


dz 

= -(d - 2) z 0 + N + 2 - 6z 0 g K , 

k—— = (d — A) g K + 2(N + 8) g 2 . (37) 


The equation for g K defines the usual one-loop /3-function; in this approximation this 
equation decouples and can be solved explicitly: 


9k 


1 + 



(38) 


where g* is the value of g at the IR fixed point, g* = (4 — d)/(2(N + 8)), and n c the value 
of k for which g K = g*/ 2. We have (g* g\): 


^c\ 

A) 


d -4 


9-9 A 
9A 


<£_ 

9 a’ 


(39) 


= uKd/(3g*) is the typical scale which separates the scaling region, dominated by 
the IR fixed point, where g = g*, from the perturbative region, dominated by the UV 
fixed point g — 0 (when k n c , one can expand g K in powers of n c /n; in leading order 
9 « = {u/ 3)(l-K/«:) 4 - rf )). 

We show in figs. 0 and m the dimensionless coupling g K and the anomalous dimension g K 
obtained by solving the complete LPA’ equations numerically for d = 3 and N — 2. The 
coupling constant g\ has been fixed to a small value, and rh\ has been adjusted in order 
to reach the IR fixed point as n —> 0. Note that w K depends a priori on u, k, and A, but 
since it is dimensionless, it can only depend on the ratios k/u and uj A (in d = 3). However, 
because the theory is super-renormalisable in d = 3, w K becomes independent of u/A in the 
limit of large A. One finds numerically that this regime is attained when u/A ^ ICR 3 , a 
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FIG. 5: The dimensionless coupling g K as a function of n/u (in a logarithmic scale) obtained 
by solving the LPA’ equations for N = 2 and d = 3. The value of g K at the IR fixed point is 
g* = 0.064. The value of k for which g K = g */2 is k c = 0.072. 



FIG. 6: The anomalous dimension g K as a function of k/u (in a logarithmic scale) obtained by 
solving the LPA’ equations for N = 2 and d = 3. The value of g at the IR fixed point is g* = 0.044. 

condition satisfied in all numerical results presented in this paper: more precisely, we used 
u /A = 3 K 3 -1 1CT 6 ~ 1.8 1CT 4 , i.e., g K= a = 1CT 6 ; the corresponding value of m\ needed to 
reach the fixed point is rh\ = —3.999527 • • • x 10 -6 . The general behaviors seen in figs. Eland 
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El are those expected from the approximate analytic solution discussed above, in particular 
the fixed point values reached at small k. On a logarithmic scale, the change of regime 
between the perturbative regime at large k, and the scaling regime at small k occurs rather 
rapidly, at the typical scale k c . In fig. 0 k c /w ~ 0.07, not far from the value obtained in 
the approximate analysis presented above: from eqs. (EH and AMD, for d = 3 and N = 2, 
k c /u = 20/i 3 /3 ~ 0.11. 

Before closing this subsection, let us write the flow equations for the 2- and 4-point 
functions in vanishing external held, in the LPA’ limit, in a form that we shall use later. 
These equations are obtained by deriving once and twice eq. m with respect to p, then 
setting p — 0, and using the definitions in eq. (OH . They read, respectively: 


k q m 2 _ (N + 2) (2) 

^ 9k, 1 d ) 


and 


= {N + S)gllf(K) - ^(JV + 4) hjf (ft), 


r( 2), 


(40) 


(41) 


where we have defined 


4”’w - 


d d q 


(2n) d 

^d-\-2—2n 


nd K R K (q 2 )G n LPA ,(K] q) 


= 2/C 


K 


zjr 1 


l - 




d + 2 


(42) 


(1 + ml) r 

the explicit form in the second line being obtained for the Litim regulator. Note that, 
after going to dimensionless variables and making the same approximations that leads us to 
eqs. (1371) (neglect the second derivative of w K with respect to z , and assume \m 2 K \ <C 1) one 
can transform eq. m into the second of eqs. (GUI). 

For further use, we also rewrite eq. m in the following form: 

Kd K g K = (N + &)g 2 K lf\K)(l-F K ). 


where 


lN + illMK = liV + 4 
“ 2JV + 8/W( K) 3 2 2N + 8 [ 


h 

2\ ' l k 


(43) 


(44) 


The function F K gives a measure of the relative magnitude of the contribution of the 6-point 
vertex term in the flow equation for the 4-point function. One can see on fig. |71 that, as 
expected, the relative contribution of the 6-point vertex is negligeable in the perturbative 
regime (k k c ), but becomes of order 1 in the scaling regime (k -C k c ) . 
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FIG. 7: The function F K in eq. (1431) as a function of k/u (in a logarithmic scale), calculated for 
N = 2 and d = 3. 

D. Correlation functions at large N 


In the critical case, the derivative expansion gives accurate results for the correlation 
functions and their derivatives only at zero external momenta. In order to get insight into 
the effect of non vanishing external momenta we consider now the correlation functions in 
the large N limit (at fixed uN). Our goal here is to illustrate some general features of 
the momentum dependence of the correlation functions, and how this is affected by the 
regulator, not to present a consistent discussion of the flow equations and their solutions, 


which can be found in the literature 111, 2 


33i. [m ] . Thus we shall not attempt to solve 


directly the NPRG equations: since they do not close, their solution requires a somewhat 


elaborate treatment (see e.g. |22l]). Rather, we shall simply write the solution for the first 


n-point functions, relying on well known results j^j], and verify that they do satisfy the 
NPRG equations. 

For vanishing held, the inverse propagator is of the form 


G (k;q) = q +m K + R K (q), 


(45) 


where m K is a running mass given by a gap equation 


m K = r + 


Nu 

IT 


d d q 

(2vr) Q 


(G(«;g)-G( A;g)). 


(46) 
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The 4-point function has the following structure: 


r i 4 234( K ;Pl>P2,P3,P4) = S 12 5 U g K (p 1 +p 2 ) + <5l3<524^«(Pl + Pi) + ^14^23^«(Pl +P 4 ), (47) 


where g K (p) is given by 


9k{v) 


u 1 

31 + ^B„( Ki p)’ 


(48) 


with 

B d (n;p) = J -^ 1 G(n-,q)G(K]p + q). (49) 

Finally we shall need shortly the 6-point function I'f 23 4 rnrn ( K: i Pi, P 2 , P 3 , P 4 , g, — 9 ) (summation 
over repeated indices is understood) 

^7 r S 6 234mm(«; Pi, P2, Pi, P4, 9, -Q) 

= h K (p 1 + P2)5i 2 Su + h K (pi + p 3 )5 13 S 24 + h K (pi + p 4 )5i4<^23, (50) 


with 

Kip) = Ng K (0)g 2 K (p) J G 2 (K;q)G(K-,q + p ). (51) 

All these results can be obtained in a straightforward fashion by calculating the corre¬ 
sponding Feynman diagrams with a regulator. It is however easy to verify that the various 
n-point functions that we have just written are indeed solutions of the flow equations in the 
large N limit. 

To this aim, one notes first that eq. © reduces to an equation for the running mass: 

d K m 2 K = ~Ng K (0) J -^-d K R K {q)G 2 {K-q), (52) 

and using eq. (HH1) . it is easy to check that this equation is equivalent to the gap equation, 

eq. m- 

Next, we observe that in the large N limit, a single channel contributes in eq. (1111) for 
the 4-point function; one can then use the following identity in this limit: 


I ' { i2 i j(^Pi,P2,q,-q-Pi-P2)^ 3 4ij(^p 3 ,P4,-q,q-p 3 -pd = Ng 2 ^ + p 2 )S 12 S 34 , (53) 


44) 


together with eq. © for r^ 6 ), and obtains: 


nd K g K {p) = Ngl{p)jf\^p) - ^h K (p)lf ] (k), 


(54) 
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where the function I^\k) is that defined in eq. (14211 . with here n — 2 and the propagator 
© replacing Glpa 1 - The function J^\k;p) is obtained from the general definition 

j{ d\ K VP) = J 7 ^^^d K R K (q)G n ~ 1 (K]q)G(K]p + q). (55) 

Note that J^\n]p = 0) = I^\k). Explicit expressions for the function J^\k;p) are given 
in app. m 

At this point we remark that the flow equation for g K (p) can also be obtained directly 
from the explicit expression (148(1 . in the form: 


d K g K (p) = ~gfc(p)d K J 


■G(k\ q)G(n-,q + p). 


(56) 


It is then straightforward to verify, using eqs. © and (EH) that eqs. © and (© are indeed 
equivalent. The first term in eq. © comes from the derivative of the cut-off function in 
the propagators in eq. ©, while the second term, which involves the 6-point vertex, comes 
from the derivative of the running mass in the propagators. 

Note that eqs. © for m K and (15411 for g K (p = 0) become identical respectively to eqs. (1401) 
and © of the LPA in the large N limit, a well know property 0. 

In view of the approximations that we shall develop in the next section, it is worth 
analyzing characteristic features of the function g K (p). For simplicity we specialize for the 
rest of this subsection to d = 3. Furthermore, for the purpose of the present, qualitative, 
discussion, one may assume m K = 0. This allows us to obtain easily g K (p) from eq. (1481) in 
the two limiting cases k = 0 and p — 0. In the first case, we have 

1 


/ \ u 

9k( 0) = o 


31 + f^' 


(57) 


This is identical to eq. ©, with here g* = 1/(2N) and n c = Nu/9ti 2 . (The corresponding 
expressions for eq. (|3%j) involve A^ + 8 instead of N, so that the values of g* and k c obtained 
in the large N limit may be numerically quite different from the actual LPA values when N 
is not too large, e.g. when N = 2). In the other case, we have 

1 up 


u 


g K =o(p) = o 


3l + ^p 3 p + Pc 


(58) 


with p c = uN/48. 

One sees on eqs. © and © that the dependence on p of g K=0 (p) is quite similar to 
the dependence on k of g K (p = 0). In particular both quantities vanish linearly as k —» 0 or 
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FIG. 8: The function g K (p ) (in units of A) obtained from a complete numerical solution of eqs. (1521) 
and (EU, as a function of k/u (in a logarithmic scale) for five values of p: from bottom to top, 
p/u = 0.001,0.01,0.1,1 and 10. The envelope corresponds to p = 0. This figure illustrates the 
decoupling of modes: for each value of p, the flow stops when k ^ ap. The various horizontal 
asymptotes (dotted Hines) correspond to the single value a = 0.54. 


p —> 0, respectively. The result of the complete (numerical) calculation, including the effect 
of the running mass ( i.e., solving the gap equation (l-EHl) and calculating g K (p) from eq. (|48|) ). 
can in fact be quite well represented (to within a few percents) for arbitrary p and n by the 
following approximate formula 


9*(p) 


u X 


X = 


K 


p 


(59) 


3 1 + X K c Pc 

This simple expression shows that p, when it is non vanishing, plays the same role as k as 
an infrared regulator. In particular, at fixed p, the flow of g K (p) stops when X becomes 
independent of k, i.e., when k ^ p(ft c /p c ), with K c jp c = 16/37T 2 « 0.54. This important 
property of decoupling of the short wavelength modes is illustrated in fig. HI As shown by this 
figure, and also by the expression El, the momentum dependence of the 4-point function 
can be obtained from its cut-off dependence at zero momentum. In fact fig. H suggests that, 


20 












to a very good approximation, there exists a parameter a such that g(n;p) ~ 0) when 

n > ap , and g(/qp) ~ g(n = otp\ 0) when k < ap. From the discussion above, one expects 
a ~ k c /Pc — 16/37T 2 as 0.54, which is indeed in agreement with the analysis in fig. |S1 



/o\ /q\ 

FIG. 9: The function (k', p)/1^ (k) as a function of p/n (in a logarithmic scale), for different 

values of k: k = 10~ 3 ri (circles), k = u (diamonds) and k = 10 3 tt (squares). 


In order to understand better the origin of this result, we rewrite eq. (151) as follows: 

d K g K (p) = N gl(p)jf ] (K]p)( 1 - F(tt,p)) (60) 


where 

= 1 K(p)ld\ K ) 
2 


(61) 


When p = 0, eq. flUD coincides with the LPA equation f£3D> and the function F[n]p) with 
the large N limit of the function defined for the LPA in eq. The p-dependence of 

J^\k;p) is relatively simple: when p -C k, J^\k;p) ~ I^\k)] when p > k, J^\k]p) 
vanishes as 1/p 2 . On a logarithmic scale the transition between these two regimes occurs 
rapidly at momentum p ~ k, as illustrated on fig. EH Fig. HI shows the similar behavior of 
the function h K (p)/h K , where h K (p) is the function (I5T1) which appears in the numerator of 
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p/K 


/q\ /q\ 

FIG. 10: The function J 3 '{k\p)/I^ (k) as a function of p/n (in a logarithmic scale) (full line). 
The function h(n;p)/h(K] 0) as a function of p/n (dashed line). 

eq. (ED- Finally, fig. fTTI displays the function F(n]p)/F K : as one can see, the momentum 
dependence of F(k;p) is non negligible only in the region where the function J^\k]p) is 
negligible, namely for k ^ p. All this suggests that one can rewrite eq. m for 9k{p) as 
follows: 

0.9. (p) « Ngl(p) 6(1 - ( K )(1 - (62) 

where a is a parameter of order unity. Eq. © is just eq. (HH1) in the large N limit, and 
for k > ap. The ©-function ensures that the flow exists only when k > ap , and stops for 
smaller values of n. These are precisely the features observed in fig. |H1 

III. TOWARDS THE SOLUTION OF THE NPRG EQUATIONS FOR ARBI¬ 
TRARY MOMENTA 

Our proposal to solve the NPRG equations for the n-point functions at arbitrary mo¬ 
menta, builds upon the lessons learnt in the specific examples discussed in the previous 
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p/K 


FIG. 11: The function F(n;p)/F K as a function of p/n (in a logarithmic scale). 

section. Namely, we shall take advantage of the decoupling of modes, exploit the solution 
of the LPA’, and use the possibility to increase accuracy through iterations. 

The decoupling of modes is well illustrated in fig.|S]of the previous section. It suggests that 
the momentum dependence of the n-point function can be deduced from their K-dependence, 
as obtained from the LPA’. To be more specific, assume for simplicity that all external 
momenta are of the same order of magnitude, and call them generically p. Then as long as 
hi > p, one can use the LPA’ to calculate the n-point functions. When hi ^ p, the flow stops 
and the n-point functions remain at their values for hi ~ p. Note that this argument ceases 
to apply when the momenta enter as exceptional configurations, for which effectively p = 0. 
These exceptional configurations cause special difficulties that we shall have to deal with. 

The possibility to increase the accuracy through iterations is based on the property 
recalled in the previous section, that the iteration of the NPRG equations, starting with 
the classical values of the n-point functions as initial input, reconstructs the usual loop 
expansion. Thus, one may expect to improve the accuracy of the n-point functions at high 
momenta by iterations. The situation at small momenta is more subtle. Indeed, in the 
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critical regime, iterations may affect the fixed point structure, and may result in unphysical 
behaviors. This particular feature will be discussed in 0. 

The procedure that we propose starts with an initial ansatz for the n-point functions 
to be used in the right hand side of the flow equations. Integrating the flow equation of a 
given n-point function gives then the leading order (LO) estimate for that n-point function. 
Inserting the leading order of the n-point functions thus obtained in the right hand side of 
the flow equations and integrating gives then the next-to-leading order (NLO) estimate of 
the n-point functions. And so on. 

The equations will be solved starting at the bottom of the hierarchy, that is, with the 
equation for the 2-point function. The flow equation for the 2-point function involves in 
its right hand side the propagator (hence the 2-point function), and the 4-point function. 
To determine the 2-point function in leading order, we need therefore an initial ansatz for 
the propagator and the 4-point function. Similarly, to get the 4-point function in leading 
order, we need an initial ansatz for the propagator, for the 4-point function and the 6-point 
function. And so on. 

There is no small parameter controling the convergence of the process, and the terminol¬ 
ogy LO, NLO, refers merely to the number of iterations involved in the calculation of the 
n-point function considered. Obviously, the calculations become increasingly complicated 
as the number of iterations increases, and it is essential that the initial ansatz be as close 
as possible to the exact solution so that only one or two iterations suffice to get an accurate 
result. Our main task then is to construct such a good initial ansatz. 

A. The construction of the initial ansatz — Generalities 

The initial ansatz for the n-point functions are the solutions of approximate flow equations 
obtained by making the following three approximations. 

1) Vertices are slowly varying functions of the external momenta 

Our first approximation («4i) exploits a crucial property of the NPRG: the derivative 
d K R K (q) limits the range of integrations in the flow equations to q ^ hi. The momentum q 
enters the vertices in the flow equations typically in the form T^ .. n (^',Pi,P2, ■ ■■,Pn-i + q,Pn — 
q). Approximation A\ consists then in assuming that, for any set of external momenta 
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{p 1 ,p 2 ,...,p n }: 


I 'l2...ni^PhP2,--,Pn-l + g, Pn ~ g) ~ „ n (/C; Pi, p 2 , • • •, Pn-1, Pn) 

r i2 ) ..n( K ;Pl>P2,--,Pn-l,Pn) 

This approximation is justified when the momenta {pi,P 2 , ■■■,Pn} are much larger than k, 
since then we can neglect q compared to Pi, assuming that T^ is a smooth function of the 
momenta when these are large. In the opposite case of vanishing p^ s, we use the fact that 
the regulator insures that T^ remains a smooth functions of its arguments; in this case, the 
approximation A\ is analogous to the leading order in the derivative expansion, i.e., to the 
LPA, known to be a good approximation. 

The approximation A\ is used to set q = 0 in the vertices T^ and to factor them out of 
the integrals in the r.h.s. of the flow equations. 

2) Propagators 

The second approximation (. 4 . 2 ) concerns the propagators in the flow equation, for which 
we make the replacements: 


< 1. (63) 


G(p + q) —> GLPA'(q ) © ( 1 


a 2 p 2 


(64) 


where a is an adjustable parameter. A motivation for this approximation is the decoupling of 
high momentum modes in the flow equations, as illustrated in sect. ITTPl and the parameter 
a will play here a role similar to that it plays in sect. Ill Dl A measure of the quality of this 
approximation is provided by fig. El which shows the ratio J^ 3 \k; p) / J^ 3 \k; p = 0) where 
J^ l \k;p) is defined in eq. (IH51 1. J^ 3 \k',p) is the integral wich remains in eq. (fTTl) after 
approximation Ai and after chosing as propagator that of the LPA’ (the consistency of this 
choice will be verified shortly). As seen in fig. El J^ 3 \k; p)/J^ 3 \k; p = 0), as a function of 
p 2 /n 2 , looks indeed like a step function, with a weak residual n dependence. 

Different criteria can be used to fix a. One may fix a so that the inflexion point of the 
curve in fig. Elis at ap = tz. One then obtains, for A^ = 2, a ~ .9. One can also adjust a 
so that the integral over k of J^ 3 \k;p) is identical to that of J^ 3 \k', 0)0(1 — Oi 2 p 2 /n 2 ). This 
yields a ~ .6. We regard these two possible choices as extremes and we adopt the value 
a = .75 ± .15 for our leading order estimate in the case N = 2. 

Before moving to the next approximation, let us write the equations for the 2-point and 
4-point functions obtained at this stage, i.e., after approximations Ai and A 2 - The equation 
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for the 2-point function becomes 

Kd K r { S(K ;p) = -ir [%(k;p, -p, 0,0 )I ( f\n) , (65) 

and that for the 4-point function reads 

^« r i234(PbP2,P3,P4) = Id\ K ) 

X |© [k 2 - a 2 (pi + P2 ) 2 ) r^-(pi,p 2 , 0, -Pi ~ P2)r^(p 3 ,P4, 0, -p 3 - Pi) 

+ 0 (n 2 - a 2 (p\ + p 3 ) 2 ) r^(pi,p 3 ,0, -Pi ~ P 3 )T { 2 %(P 2 ,Pi, 0, -p 2 ~Pi) 

+ © (k 2 - a 2 (pi + p 4 ) 2 ) r^(pi,p 4 ,o, -pi -P4)r^(p 3 ,P2,o,-p 3 -P 2 )} 
-^ r f234n(Pl 5 P2, p 3 , P4, 0, 0) /^(k), (66) 

where the function I^\k) is dehned in eq. 

Note that the approximation A 2 amounts to truncate severely the high momentum tails 
of the propagators. This will cause inaccuracy at high momenta, and a dependence of the 
leading order results on the value of a. 

3) Approximation for the (n+2)-point function 

In order to close the equation for T (ri ) we need an approximation for r( n+2 \ Namely, we 
need an approximation for T (4 ^ in the equation for T^ and for ffi 6 ) in the equation for ffi 4 -*. 
Note that we do not want to perform a truncation of the hierarchy, as often done, by setting 
to zero the higher order n-point functions: indeed in the scaling regime the contributions 
of all vertices are of the same order of magnitude. Rather, we shall try to obtain a rough 
estimate for T (ri+2 ), which is sufficient in order to get the initial ansatz for (this rough 
estimate is not to be confused with the initial ansatz for ffi n+2 )). 

In order to construct this estimate for ffi n+2 ) we rely on the LPA’ and also use an ap¬ 
proximation inspired by the analysis of the correlation functions in the N —> 00 limit of the 
previous section (sect. Ill Dl) . We consider explicitly here the equations for the 2-point and 
the 4-point functions. 

In the case of the equation for T (2 ), one needs an approximation for T^(«;p, —p, 0,0), 
as can be seen in eq. m- This can be hinted from eq. (BED, leading us to assume 
Ti 2 «( K ;P, — Pi 0, 0) = N5i2g K {Q). The resulting initial ansatz for T (2 ^ is simply a momen¬ 
tum independent function, the running mass, whose flow equation is given by eq. U52D- 
Therefore, our initial ansatz for the propagator is consistent with eq. (12411 . to within the 
small effect of the anomalous dimension which is ignored at this stage. 
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We turn now to r (4 h As can be seen in sect. ETUI after doing approximation A 2 , the 
two types of terms in the right hand side of the flow equation of Id 4 ) are proportional, 
with a coefficient that depends only on k (see eqs. (HuH) and (E2), and fig. fTTl) . Onr third 
approximation (*4. 3 ) consists in assuming that this property holds in general, i.e., we set: 

3* r m4 (Pi>P2> P3,Pa) = -F K d K T$$ (pi,p 2 , P 3 ,Pa), (67) 

where in the l.h.s. r ^34 is the 6-point vertex contribution to the flow of T 1 - 4 ) (last line in eq. 
m and fig. H, while the term multiplying — F K in the r.h.s is that including only 4-point 
vertices (the first three lines in (16(411 and fig. |H). This relation becomes trivial in the LPA, 
i.e., when all external momenta are zero. This allows us to fix F K from eq. (®. 

Combining all approximations, one gets the following equation that needs to be solved 
in order to get the initial ansatz for 

Kd K T ( $ u (K-,p 1 ,p 2 ,p 3 ,p i ) = I^\k) (1 - F k ) 

X |© (k 2 - a 2 (p 1 + P 2 ) 2 ) T ( $ij(Pl,P2, 0, -pi - P2)r%(p3,P4, 0, -p 3 - Pi) 

+ 0 (k 2 - a 2 (pi + p 3 ) 2 ) r ( S ij (p 1 ,p 3 , 0, -Pi - P 3 )^2Aij(P2,PA, 0, -p 2 - Pi) 

+ 0 (k 2 - a 2 (pi + Pi) 2 ) rJ2j(pi,P4, 0, -Pi - Pi)r ( 3 2 ) ij (p 3 ,p 2 , 0 , -p :i -P 2 )} ■ 

(68) 

In the rest of this section, we construct the solution of this equation in terms of the 
solution of the LPA’. As a simple illustration of the method to be used, consider first the 
totally symmetric configuration of momenta: (pi +P 2) 2 = (p 1 + P 3) 2 = (pi + Pi) 2 = P 2 (and 
Pi = P 2 = p\ = p\ = 3p 2 /4). One then distinguishes in eq. (IfiHl) two regions, according to 
the value of k relative to ap. When k > ap , all the terms in eq. (jHHj) are non-zero. One can 
then verify that the LPA’ expression of r^ 4 \ i.e., that given in eq. K25D, is a solution of the 
equation. Since the initial condition at k = A has the form of eq. 425), and since eq. (EH) is 
a first order differential equation in k, the LPA’ solution is the unique solution for k > ap. 
When k < ap, the r.h.s. of eq. EH) vanishes and the flow stops. In this region, the solution 
remains the LPA’ solution, but taken at the fixed value k = ap. These are the features that 
we uncovered when we analyzed the correlation functions in the large N limit in sect. Ill Dl 

A similar separation into different regions, according to the value of k, can be done for 
general momentum configurations. In all cases, when k is larger than all the combinations of 
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momenta appearing in the ©-functions in eq. (168(1 . the solution is simply the LPA’ solution. 
The other regions, where some of the ©-functions vanish, have to be analyzed case by case. 
One can then solve eq. (EHD in two steps: first, for one vanishing momentum, p 3 = 0 and 
p 4 = —p\ — p 2l then for any combination. 

In the next two subsections, in order to simplify the notation, and except when ambigu¬ 
ities may arise, we shall often omit to indicate the explicit ft dependence of T^ 4 b 


B. 


Calculation of T 


(4) 

1234 


(Pl,P2,0, ~Pl -p 2 ) 


In this case, eq. (IfiHI) reads: 

^K r i234(Pl>P2, 0, -Pi - p 2 ) = 1 - F k ) 

x I© (ft 2 - a 2 (pi + p 2 ) 2 ) ^[%(pi,P2,0,~Pi -P2)r^-(0,-pi -P2,0,P1 +P 2 ) 

+ 0 (hi 2 - a 2 p\ ) 0, 0 , ^Pi)T^ ij (p2, ~P\ ~ P2, 0,pi) 

+ © (ft 2 - a 2 p\) r^.^pi, -pi - P 2 , 0,p 2 )r^(0,p 2 , 0 , -p 2 )} . (69) 

Notice that in each term in the r.h.s. there is one vertex evaluated with two vanishing mo¬ 
menta. Furthermore, because of the theta functions, each term gives a non-zero contribution 
only when the remaining non vanishing momentum is smaller than n/a. We are therefore 
in the conditions discussed at the end of the last sub-section: the 4-point functions with 
two vanishing momenta are simply the LPA’ ones (eq. (12511 ). By using the fact that bosonic 
vertex functions are completely symmetric under simultaneous exchange of internal indices 
and momenta, we can rewrite eq. © in the following way: 


x j© (ft 2 - a 2 (pi +p 2 ) 2 ) ^nn{pi,P 2 , 0 ,-pi -P2)^34 

+ r i234(Pl>P2,0,-pi -p 2 ) +rS? 43 (Pl>P2,0,-pi -p 2 ) 

+0 (ft 2 - a 2 p 2 ) rS 4 (pi,p 2 , 0 , ~P\ -p 2 )S 13 

+ r 3214(Pl,P2,0, -pi -P 2 ) +rS? 34 (Pl 5 P2,0, ~Pl -P 2 ) 

+0 (ft 2 — ct 2 '"' 2 '' 


l pl) r iS 4 (Pl,P2, 0 , -Pi- p 2 )^23 
+ r i324(Pl ? P2,0, -P 1 -P 2 ) + rg>4(pi, ft ,0, -Pi -P 2 ) 


( 70 ) 








This is a first order linear equation where the momenta are parameters. To solve it, we can 
assume without loss of generality that p\ > p\ > (pi+Pi) 2 - For the rest of this section, except 
when that would lead to confusion, we will drop the arguments of ^ 34 , being understood 
that T 4234 refers to T[ 4 2 ] i 4 (Kyp[,p2, 0 , —p\ — p 2 ). We need to consider four different regions, 
according to the value of k: 

a) k > ap±. In this region, the solution is identical to that of the LPA’. 

b) api > k> ap 2 - In this region, eq. (0 becomes: 

K 9 «r < 4 2 > 34 = g jf\n)( 1 - f k ) {rgpsi + 2rf 2 ) 34 + rfi + r[%s 23 + rg*}. (71) 

To solve this equation, we first notice that the solution is symmetric under the exchange of 
the 2nd and the 4th internal indices (with no exchange of the momenta), i.e.: 

r 1234(Pi? P2j 0 , p\ -P2) = T[% 2 (p 1 ,p 2 , 0 , -pi -p 2 ). ( 72 ) 


This property is true for k = api, and one can verify that it is maintained along the flow. 
We then look for the general solution symmetric in the indices 2 and 4, in the form: 

Fl234 = (^12^34 + hl4^32)FA + 5l3^24Fs. (73) 


Substituing (d in eq. d, one finds the following system of linear equations: 

I Kd K r A = lf( K )g K (l-F K )((N + 4)r A + 2r B ) 

\ Kd K r B = lf(K)gJl - F k ) (2T a + 2T b ) . 

The matrix 


( AT + 4 2 

y 2 2 


has the eigenvalues (which are both positives for N > —2) 

N + 6 ± VN 2 + 4IV + 20 


corresponding to the following eigenvectors: 



(74) 


(75) 


(76) 


(77) 
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Using these eigen-vectors, one can write the general solution of eq. m as: 



where a± verify: 


nd K a^ = lf\n)g n ( 1 - F K )\±a± 


A± 

iV + 8 


Kd K (\og g K )a±. 


(79) 


We used eq. dED to obtain this result. The equation above can be integrated analytically, 
to give: 


a K — a api 


9api 


A ± 
N +8 


By imposing continuity between the two regions (a) and (b) at k 
solution in the region api > k > ap 2 : 


(80) 


api, we obtain then the 


p(4) _ 

1 1234 ~ 


9api 


A_ - A, 


<5 1.3 <5 


+ (<5l2<5 34 + ^14^23) 


13«24 


-A, 


(A- — 4) 


(—) 

\ 9api J 


9ctpi 


N+8 


+ A_ 


A + 

N+8 


9k 

9ap\ 


(A+ - 4) 


9k 

9ap\ 


N+8 


X_ 

iv+8 


(81) 


c) ap 2 > k > a\pi + P 2 1 - hi this region, eq. (17(11) becomes: 

s«r< 4 2 >34 = jU {rg/34 + r< 4 2 > S4 + rfi} a. (logs,). (82) 

We need now the general tensor decomposition: 


Ul234 — A + (<5l3<524 + ^14(532)1^5 + (<5l3<^24 — <5l4<532)rc. (83) 


By substituting in eq. (jH2j) we get: 

' d K r A = J + 8 ((N + 2)r A + 2r B )d K (\ogg K ) 

< d K T B = j^T B d K {\ogg K ) (84) 

d K r c = 0. 

The antisymmetric sector (Tp) is decoupled. In order to get the solution in the symmetric 
sector (T^jTb), we diagonalize the matrix: 


/ N + 2 

V 0 



(85) 
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and get the eigenvalues 


/I- 1_ — N + 2 , /i_ — 2 


( 86 ) 


corresponding to the eigenvectors: 

dd 1 

0 j \-N /2 

One can then write the general solution of the symmetric part of eq. 


(87) 


as 


£ 1 


N/2 


( 88 ) 


where bt verihes 


which, using eq. (TIHl) . leads 


adiib* = Id\a)g K (l - F K ) n± 6*, 


to 


ac ap2 » » 


\ n / 

\ y&p2 / 

Imposing continuity between the regions (b) and (c), at ft = ctp 2 , one obtains finally: 


(89) 


(90) 




r (4) e /,+ A 9k \ N+8 , ( 9k \ N+8 

M234 — <^12^34 I b ap2 - j +b ap2 


9ap 2 , 


7V / n \ N +8 

- (<fl3^24 + ^14^32) ~) 

^ \ 9 otp 2 J 

+ (<^13^24 — ^14<^32)rap 2 , 


(91) 


where 


7 + _ 9api 

ap2 A_ - A h 


d A+ + >- A+ - 4) )(t 


*+ 
N +8 


b a p 2 


•(-A- 4 (W-- 4) ) (^) 

(A.-A+-4) 

\9api / 


A_ 1 
iV+8 


9ap\ 


JV+8 


IV A_ - A 

— (A + — A_ — 4) 

V fl^pi , 


A_ ‘ 

77+8 


(92) 
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( 93 ) 


'h,„ \ ”+' _ »+* 

9api / \ 9ap i / 

d) «|pi +P 2 I > K- In this region the flow simply stops. The result is then: 

r i234(«;Pi ) P2,0,-pi - Pi ) = T^fa = «|Pi +P2|;Pi,P2,0,-pi -p 2 ) (94) 

C. Calculation of — p, q, —q) 

At this point, we could solve eq. (IfiHj) for any combination of momenta, given the fact 
that once the function is known for the particular combination that has been treated in 
the last section, all the information appearing in the right-hand side of the equation is 
known. To give an example we shall consider in this sub-section the explicit calculation of 
Fit —PiQi ~Q)i when q < k. The result will be used in the next section as the initial 
ansatz for T^ in the calculation of the LO expression of the self-energy. Note that since the 
favored values of a are smaller than 1, aq < k. 

For the considered values of momenta, eq. (EE} becomes: 

l 2 h(p> -P, <?, -q) = 4 3) (k)(1 - F k ) {r^.(p, -p, 0, 0)T^.(g, -q, 0, 0) 

+0 j ^[%(p,q,0,-p - q)T^ ) ij (-p,-q,0,p + q) 

+0 (l -a 2 <<P ^ "j T^.(p,-g,0,-p + g)r^.(-p,g,0,p-g)|. (95) 

The r.li.s. of this equation includes the expressions of T^ that have been determined in the 
previous sub-section. It is useful to separate the contribution to T^zqp, —p, q, —q) coming 
from the first line from those coming from the second and third lines of eq. El above. The 
first contribution corresponds to the s-channel, the second corresponds to the sum of the t 
and u channels (see fig. m 

1) s-channel 

In the s-channel, there are two kinematical regions: 

a) ap < k. In this case, we have 

K9„rW„(p, - p , q , -q) = /f( K )(l - F K )gl(N + 2)=<5 12i (96) 

whose solution is (see eq. dH): 

r i2iZ (p> “P> -q) = 9j 12- (97) 


pC _ 9ap\ N + 2 

" P2 ““ A_ — A+ 2 
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b) ap > k. In this region, we have: 


-P> -q) = 4 3) ( k )(! - F *)( N + 2) 2 5l2^ap ( — 

\9ap 


IV+2 

IV+8 


(98) 


whose solution is: 

r i2w(P> 9> -9) = 9ap^l2 (A + 2) 


A + 2 / 

A + 8 V^p 


JV+2 

IV+8 


(99) 


where we used continuity between the two regions in order to fix the integration constant. 

2) t and u channels 

Let us now turn to the contribution of the t and u channels in eq. ©. Since the two 
channels only differ in the sign of q, we consider only the ^-channel, with kinematical variable 
\p + q\. There are two situations to analyze. 

A) \p + q\ > \p\. In this case, there are two kinematical regions: 
a) a\p + q\ < k. In this region the contribution in eq. (1(tH1) reads: 


«0* r iL(P> -P. 9, ~q) = Id’WO- ~ F K )gl 3(N + 2)5i2- 
The solution is easily obtained by using eq. d3SD in order to eliminate F K . One gets 


_ t ( 3 ) 


( 100 ) 


rfita -p.9, - 9 ) = 3 ^, + g g^i 2 - 

b) a\p + q\ > k. In this region the flow stops and we obtain: 


N + 2 


( 101 ) 


r Su(P. -P.9, -9) = 3 ^, + g g»lp+a|fe 


iV + 2 


( 102 ) 


B) \p + q\ < p. In this case there are three kinematical regions: 

a) ap < k. This region is identical to (Aa) above. The solution is given by eq. (jlOljh 

b) ap > k > a\p + q |. Here the contribution to eq. (1HH1) becomes : 


«^« r t lu(p, -p, q, -q) = F d\ K )( l ~ F *)5i2 


— r( 3 )( 


N + 2 2 

N 2 + m + 20 9ap 



-N + 6N + 30 + 


N+ 14 


VN 2 + AN + 20 ) ( — 

2 /V (Jap 


2A_|_ 

lV+8 


+ -N 2 + 6N + 30 


N+ 14 


VN 2 + AN + 20 ) ( — 
2 / \ 9 exp 


2A_ 

77+8 


( 103 ) 
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and has as solution: 


r iL(P> “P> -9) = + <$12 


V + 2 


iV 2 + 4iV + 20 


Sap 


X 


fiV 2 + 6V + 30 + A±iVjV 2 + 4V + 20 


2A . - N - 


9k 

9ap 


2X + -N-8 
JV+8 


- 1 


| N 2 + 6N + 30 - ^hyFT4f+20 
2A_ — N — 8 


9ap 


N+8 


- 1 


(104) 


where, again, we have imposed continuity. 

c) k < a\p+q\. In this regime the flow stops. The solution is found by fixing k = a\p + q\ 
in eq. (110411 . 

IV. LEADING ORDER RESULTS 


The self-energy £(fc;p) is obtained by integrating eq. 0 from the microscopic scale A to 
the given value of the parameter k: 

5i 2 Z(k;p) = 5 12 r - ^ dn' J ^^G 2 (k/; q)d Kf R K >{q)T^ ii (K , -,p f ^p, q, -q), (105) 

where we have used the boundary condition E(/c = A;p) = r, with r the bare mass. We shall 
be working in the critical regime, i.e. for a vanishing physical mass. Thus r is supposed to 
be adjusted so that T,(n = 0;p = 0) = 0, that is: 

l fV t diq 


^12 r = -r / dh J I T ^-G 2 (K'- 1 q)d K/ R K '(q)T ( $ ii (n , - 1 0,0,q,-q). 


/0 J (27T) 


(106) 


One may use this equation to eliminate the explicit r-dependence in eq. m 

fa S(k;p) = J o (1k ' j Mifa'lZh-Zhe,-?) 

+ \f Q dK> j j^ G2 ( rd 'Gi) d K' R K'(q) -q) -rSU^; 0 , 0 ,?,-?)), 

(107) 

from which one immediately deduces the following expression for the physical self-energy 
£(p) = E(/c = 0;p): 

5i 2 E(p) = ^ [ t^ = ^G 2 (k / ; q)d K >R K/ (q) (t^k'^p, -p, q, -q) - rj?*(«'; 0, 0, q, -q)') . 


(2vr) c 


(108) 
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This expression automatically satisfies the criticality condition at £(p = 0) = 0. But, of 
course, it holds provided eq. m holds. 

In the following subsections we study the self-energy at leading order (LO) of our ap¬ 
proximation scheme. The leading order consists in using in the r.li.s. of eq. cnzD the initial 
ansatz for T^tX; P, —p, q, —q) that has been derived in the previous section. Note that for 
this initial ansatz, T^fV; 0, 0, q, —q) is given by the LPA’ expression, so that eq. (110611 is 
satisfied at LO with the value of r obtained by solving the LPA’ (eq. (110611 for r is then 
equivalent to eq. PH). a self-consistent equation for the running mass m K where the value 
of r is adjusted so that, for a given value of the bare coupling, rn K = 0)). 


A. The self-energy at LO 


As we just mentioned, in order to calculate £lo> we use, as input in the r.h.s. of eq. PHD, 
the initial ansatz for both the propagator and the 4-point function. The initial ansatz for the 
propagator is needed only for q < k and is taken to be the LPA’ propagator (see eq. (J23D) : 

G lp A <(«; q < «) = Z K q 2 + ml + R K (q ) = Z k k 2 ( 1 + mjj). (109) 

The initial ansatz for T^ was determined in section imn It depends only on k, p 2 , q 2 and 
the angle 6 between p and q. By performing the integrations over all angles other than 9 , 
one gets: 


d - 1 r K 

<5i 2 £lo(«;p) = — ^ — K(l ~ l I dK ' 


q d 1 dq I 2 + q K i{^- - 1) 


o ““ Z k ,k' 2 (1 + rh 2 K ,) 2 

x [ d9sm9{l-cos 2 e) {d ~ 3)/2 T ( ^ ii (K , ;p,-p : q,-q) 

/ 

1 


' K 


/2 


LAav, f d K ' 


Z k >k ,2 ( 1 + m 2 ,) 2 


q dq[ 2 + Tj K >(—-l) 


K 


x / d6 sin 6(1- cos 2 6) (d 3)/2 -p,q, ~l) ** r(h(«';0,0, q, -q) 


( 110 ) 


The physical self-energy at LO is then given by: 


&12 £ lo ( p ) — — Kd~ i 

47T 


dn- 


q d ~ l dq[ 2 + rj^-l) 


Z K n 3 (l + m 2 K ) 2 J 0 

x f d6 sin 6(1 - cos 2 9) (d ~ 3)/2 (V HU^p, -p,q, -q) - r^/c; 0, 0, q, -g)) . 


(Ill) 
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This expression has interesting scaling properties that we shall present for the case d = 
3 (most of the discussion extends to arbitrary dimensions, with the replacement of u by 
M V( 4 —d))_ 

First, a simple analysis shows that T,lq(p) in eq. (11111) can be written in the form 
E lo(p ) = w 2 E(p/w) where E is a dimensionless function. To see that, we note that, as 
seen in section UlI Cl F^ui^P, — P, q, — (?) is proportional to the LPA’ function gi where 
l = k, ap, or a\p + q\. Now, as discussed in section 111 Cl in d = 3, the dimension¬ 
less function gi l d ~ A gi only depends on l/u if u/A is small enough. It follows that 
—p, q, — q ) = uF^u(k/u;p/u, —p/u, q/u, — q/u ) where f ^ is a dimensionless func¬ 
tion. The result for E^o follows after noticing that the remaining dependence in A sits in the 
upper limit of integration: since the integral converges, that dependence becomes negligible 
when A/« > 1. 

A similar (but approximate) scaling holds for the dependence on the parameter a. To see 
that, let us set q = 0 in the 4-point functions in eq. (11111) (similarly to what is done for the ap¬ 
proximation A\ of sect ion ITl I Al) . Then, by using the explicit expressions of T (4 ) (/t; p, —p, 0, 0) 
nresented in sub-section Ilf f Cl one can verify from ea. m that E lo(p) is a function of 
ap only, i.e., E lo(p) = E (ap). In fact, we expect this property to be best satisfied for low 
values of a: Indeed, since the second line of eq. cm is non vanishing only for q < n < ap 
(see sub-section unni). the smaller the value of a the smaller the domain of variation of q, 
and the better is the approximation q = 0. The approximate scaling on a is clearly visible 
in fig. [HI 

Turning now to the momentum behavior of E lq, we note that both the low and high 
momentum regimes are correctly reproduced, independently of the value of a. At large mo¬ 
menta, we recover the logarithmic behavior predicted by second order perturbation theory, 
namely E (p) ~ In (p/u). However, the numerical coefficient in front of the logarithm (which 
does not depend on a) comes about 7% higher than the correct one (( N + 2)u 2 /{288ir 2 )). 

In the low momentum region, we obtain the expected power law behavior p 2 + E lo(p) ck 
p 2 ~ v * . It turns out that the exponent g* is the value of the function g K at the IR fixed point of 
the LPA 7 . This is verified numerically with a numerical uncertainty of 0.001, independently 
of the value of the parameter a. But this is also an exact result at the present level of 
approximation. To prove this, let us first note that, in eq. cm, the difference of the 
two functions H 4 ) in the second line is non-vanishing only when n ~ p. Indeed, as can 


36 










FIG. 12: X(p) (in units of A 2 ) as a function of otp/u for N = 2 and various values of a: a = 0.6 
(circles), a = 0.7 (square), a = 0.75 (diamond), a = 0.8 (triangle up) and a = 0.9 (triangle left). 
The curves exhibit the a scaling explained in the text. 


be easily seen from their explicit expressions given in section mn the two functions R 4 ) 
coincide when k > ap and k > a\p ± q\] therefore there are two different situations where 
the contributions are not zero. The first situation is k < ap. This implies that k < p 
(remember that a < 1). The second situation is more subtle. If we have k < a\p ± q |, one 
has \p\ > \p ± q | — |g| > \p ± q\ — K > n(l/a — 1), where we used the triangular inequality 
and the fact that q < hc. In both situations we found that, as announced, the integrand in 
eq. (jnu> is non vanishing only when k < (3p , where (5 is a number of order 1. It follows that 
if p is in the scaling region, i.e., if p -C p c (with p c ~ k c ), so are all the momentum variables 
in the integrand of eq. m, i-e-, p, q, \p ± q\, and tz. Then all the functions appearing in 
the r.h.s. of eq. «nu> are in the scaling regime, and their dependence on k is controlled by 
the IR fixed point: 


ml ~ rh* 2 , g K ~ g* , r) K ~ rj* , Z K (x k v * , 
where we used eq. m- Then, from eq. flu , we get 

9 » CC 


( 112 ) 


(113) 
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At this point we perform the change of variables k = px and q = py in order to make explicit 
the p dependence of E(p) in eq. (II 111) : we collect a factor p v * from Z~ x and an overall 
factor p d ~ 2 due to the terms dn, dq, «C 3 and q d ~ 1 appearing in the integrand. As for the p 
dependence of the 4-point functions, one uses the fact that T^ 4 ) («; p, — p, q, —q ) is proportional 
to gif(g K /gi) where l is either ap or a\p ± q |, and / a dimensionless function (see section 
urn After the change of variables, using eq. m , r(4) («; p, -p, q, -q) - r(4) («; p, -p, o, o) 
can thus be written as p A - d ~ 2r i times a function of x and y. Altogether, and using the fact 
that, as shown above, T, LO (p) — w 2 E(p/w), one gets: 

Z LO (p)=Cu^p 2 - r >\ (114) 


where the proportionality coefficient C is the remaining dimensionless and finite integral over 
x, y (and 6), which only depends on the parameter a. 

The anomalous dimension obtained from the present calculation is then identical to that 
calculated in the LPA’ (see eqs. m and (1221) 1. Its value, g* ~ 0.044 is to be compared 


with the best estimates available in the literature, e.g. g = 0.0354 ± 0.0025 


36]. A sim¬ 


ple proof that the dependence of the field renormalisation factor on the scale k determines 
in general the momentum dependence of the self-energy (thus defining the anomalous di¬ 
mension) is presented in app. El However there is no guarantee that this property should 
hold in any approximation (for instance it does not hold in the derivative expansion). It 
is therefore gratifying to see that that the power law behavior expected for the momentum 
dependence of the self-energy in the scaling regime comes out naturally in the LO of the 
present approximation scheme. 


B. Calculation of A(p 2 ) 


As a further test of the quality of the leading order result for the self-energy, we have used 
to calculate the shift A T c of the transition temperature of a dilute, weakly interacting, 
Bose gas. It has been shown that A T c is linear in an 1 / 3 Q, where a is the scattering length 
and n the particle density: 


A T c 


7^0 

^ C 


c an 


1/3 


(115) 


Here T c ° is the condensation temperature of the ideal gas and AX) 
transition temperature of the interacting system. As shown in Ref. 


T c — X c ° with T c the 
, the coefficient c can 
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be related to the change A (ip 2 ) in the magnitude of the fluctuations of the field described 
by the action ©: 

2567T 3 A {ip 2 ) 


c = 


(116) 


(C(3/2)) 4/s Nu 

in the limit u —> 0 (and for N = 2). 

The best numerical estimates for A(ip 2 ), and hence for c, are those which have been 
obtained using the lattice technique by two groups, with the results: c = 1.32 ± 0.02 


• 05 Q. 


and c = 1.29 ± 0.05 |19||. The availabilty of these results has turned the calculation of c into 
a testing ground for other non perturbative methods: expansion in 1/N jhil 3 ? 1 ] , optimized 


perturbation theory |38J, [39|, resummed perturbative calculations to high loop orders 
Note that while the latter methods yield critical exponents with several significant digits, 
they predict c with only a 10% accuracy. This illustrates the difficulty of getting an accurate 
determination of c using (semi) analytical techniques. 

To understand better the origin of the difficulty, let us write A (tp 2 ) as the following 
integral 


A(</b 2 ) = f d 3 P 

N J (2t r) 3 \p 2 + Z(p) p 2 


\ 1 

f dp 

) 2 tt 2 J 

P 


p- 


pO 


p 2 + S (p)_ 


(117) 


where £(p) is the self-energy at criticality, i.e., £(0) = 0. In eq. (11171) . the term within the 
square brakets, to be referred below as the integrand, is, to a very good approximation, equal 
to S(p)/p (one finds numerically that this is a good approximation as soon as p/u ^ 10 -5 ). 
As we shall see shortly, T,(p)/p is peaked in the region of intermediate momenta between the 
critical region and the high momentum perturbative region (see fig. ITil . The difficulty in 
getting a precise evaluation of the integral (HHD is that it requires an accurate determination 
of £(p) in a large region of momenta including the crossover region between two different 


physical regimes Q, G] • In that sense, the calculation of c can be viewed as a very stringent 
test of the approximation scheme. 

A plot of the integrand of eq. (11171) (divided by a ) is shown in fig. [T31 for various values 
of cc. As announced, the momentum at which the integrand reaches its maximum lies 
in the intermediate momentum region: in fig. EH this is ap/u ~ 0.2. The approximate 
scaling behaviour that can be observed in fig. 1 131 follows from the property of the self-energy 
discussed in sect. CO: as we have seen there, £lo(p; a ) — £(op), so that, setting p = ap 


A(^) 


N 


—a 


d 3 p £(p) 


(27t) 3 p A 


(118) 
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FIG. 13: The integrand of eq. (11171) (divided by a , and in units of A) as a function of ap/u for 
various values of a: a = 0.6 (circles), a = 0.7 (square), a = 0.75 (diamond), a = 0.8 (triangle 
up) and a = 0.9 (triangle left) (points shown are those needed to the numerical calculation of the 
integral in eq. dZJ.) The curves exhibit the approximate a scaling explained in the text. 

Fig. El shows the value of the coefficient c as a function of a. The (almost) linear 
behavior of c as a function of a follows directly from eq. (EH1). Deviations from the non 
linear behaviour can be seen for a ^ .7: as we have discussed in the previous subsection, for 
these large values of a, the approximation S LO (p; a) ~ E(ap) becomes less accurate. As we 
can see, when a = 0.75 ± 0.15, one gets c = 1.3 ± 0.3. This result confirms the quality of 
the leading order expression of the self-energy for all values of the momentum. 

We have also calculated for different values of N, and compared the corresponding 
results with those obtained by different means and available in the literature. The quality 
of our numerical estimates remains of the same level as long as N ^50, but for larger values 
of N, the calculations lose accuracy. The range of acceptable values of a (see sect. 
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FIG. 14: The coefficient c calculated in L.O. as a function of the parameter a. 


remains the interval ~ 0.6 — 0.9, and the resulting error bars on the predicted value of c stay 
of the order of 23 — 29%. One gets, for N = 1, c = 1.06 ± 0.27; for N = 3, c = 1.47 ± 0.39; 
for N = 4, c = 1.66 ± 0.44; for N = 10, c = 2.33 ± 0.60; for N = 40, c = 2.97 ± 0.63. These 
numbers are to be compared with those obtained using other methods; lattice calculation 
jdll ] or ressumed perturbation theory carried up to 7-loop order 0 give: for N = 1, 
c = 1.09 ± 0.09 (lattice) and c = 1.07 ± 0.10 (7-loops); for N = 3, c = 1.43 ± 0.11 (7-loops); 


for N = 4, c = 1.60 ± 0 
N —> oo is also known 


TO (lattice) and c = 1.54 ± 0.11 (7-loops). The exact result for 
17| : c = 2.33. One observes that, for all values of N, the best 
accepted results always lie within the error bars of our LO prediction; they approach the 
lower limit of the band when N grows (the origin of the latter property can in fact be 
understood by analyzing the steps leading to eq. 

Before finishing this section we present a consistency check of the approximation A\ made 
in sect. Erato construct the initial ansatz for the 4-point function T^b This approximation 
consists in neglecting the internal momentum dependence in the 4-point vertices appearing 
in the r.h.s. of the flow equation. This was done in order to obtain eq. (jnSl) for T*' 4 ' 1 . Here 
we shall make this approximation A\ in the equation for the self-energy, eq. fill 111 . Fig. IT51 
compares the self-energies obtained form eq. cm with and without approximation A \. One 
can see that the approximate result differs very little from the exact one. It turns out that 
both the perturbative regime and the exponent in the scaling regime are almost unchanged, 
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FIG. 15: Calculation of the self-energy (in units of A 2 ) used to test approximation A\, as explained 
in the text: the complete expression of £(p) (triangles) and the approximate one (squares). 

most of the difference being concentrated in the intermediate momentum region. This is 
verified by calculating the coefficient c with and without the approximation A \: The value 
obtained with Ai is about 10% smaller than that obtained with Elo This illustrates the 
large sensitivity of the coefficient c to variations of the self-energy in the cross-over region. 


V. CONCLUSIONS 

The calculation of the self-energy of the O (N) model demonstrates that the approxima¬ 
tion scheme that we have presented fulfills its goal, that is, it offers a simple way to calculate 
the full momentum dependence of the n-point function. The accuracy achieved in the lead¬ 
ing order is already satisfactory, over the full momentum range, as shown by the various 
tests that we performed. 

A crucial ingredient in the calculation is the construction of the initial ansatz for the 
4-point function. That in itself is an important part of the present paper. This initial 
ansatz is obtained by solving an approximate flow equation derived using well motivated 
approximations. The resulting 4-point function, albeit approximate, exhibits a realistic 
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momentum dependence, also in the entire momentum range. In particular, the power law 
behavior expected in the scaling regime is reproduced. 

The approximations that we have introduced to construct the initial ansatz for the 4-point 
function involve a parameter a that needs to be adjusted in such a way that approximate 
expressions match best the exact expressions that they are supposed to represent. This 
introduces a theoretical uncertainty, which, in the case of the calculation of the shift of the 
Bose-Eisntein transition temperature that we have presented, is of the order of 25%. 

In a forthcoming paper 2CJ, we shall present results of a next-to-leading order analysis for 
the self-energy. To do so we shall need to improve the accuracy of the 4-point function, as 
compared to the initial ansatz presented in this paper. That is, we shall calculate the 4-point 
function at leading order, i.e., construct an initial ansatz for the 6-point function. The next- 
to leading order calculation of the self-energy will allow us to test fully the approximation 
scheme, and detect some of its weaknesses. As we shall see, the calculation of c will be 
greatly improved, in particular the dependence on the parameter a will be eliminated, and 
results obtained in quite good agreement with lattice data. 


APPENDIX A: THE FUNCTION rj K AND THE ANOMALOUS DIMENSION 

It is usually accepted [fj that the ^-dependence of the held renormalisation factor Z K (for 
values of k where the couplings have approached the infrared fixed point) determines the 
anomalous dimension of the held. However it is not a priori obvious that the momentum 
dependence of correlation functions obtained after some approximation follows automatically 
the corresponding scaling law: this is not so for instance in the derivative expansion. We 
find it therefore useful to present in this appendix a simple derivation of this property. The 
arguments in the proof given here help to clarify the conditions under which this property will 
be satished in an approximation. This also completes the derivation presented in Sect. II V XI 
that the anomalous dimension calculated from the momentum dependence of the 4-function 
obtained as solution of the approximate equation derived in Sect. IHIl is indeed equal to that 
deduced from the K-dependence of Z K calculated in the LPA’. 

Let us consider the 2-point function T ( - 2 ' 1 (p, n,u) for p, k n c ~ j n order to be 
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in the scaling regime. Then, scale invariance implies that : 



(Al) 


where / is a dimensionless function of its arguments. It follows that 


r< 2 > (p, k, m) = r< 2 > (o. k, <•)/(!) 


(A2) 


where we have set f(p/n) = f(0,p/n). Note that the dependence on the microscopic pa¬ 


rameter u is entirely contained in the factor T (2 )(0, k, u) (T (2 )(0, k, u) is well defined thanks 
to the IR regulator), and the momentum dependence factors out in the scaling function 
f(p/n). This function has a Taylor expansion at small p/ k, and /(0) = 1. At this point we 
may use scale invariance again, together with dimensional analysis, in order to show that in 
the regime k -C k c : 



(A3) 


where p is constant. It then follows that for the function T (2 ^(p, k, u ) in eq. (1A2D to have a 
limit when k —> 0, we must have, for large values of p/n: 



(A4) 


where p is the same constant as in eq. (lA3j) . Thus, 



(AS) 


We can write: 


r ( 2 ) (p, K ,u) - r ( 2 ) (o, K , u) = z.p 2 + o( P 4 ). 


(A6) 


and from eq. 


r (2) (p, k,u) - r (2) (o, k,u) = (/(p/k) - i)r (2) (o, k,u). 


(A7) 


In the regime k — > 0, so that eq. H) is valid, and p ■ C k so that f(p/n) ~ 1 + C ( p/n ) 2 , 
with C a numerical constant, one can then use eqs. (I A Ml) and (na to deduce the behaviour 


oiZ K . 



(AS) 
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APPENDIX B: THE FUNCTION jf } (Ac;p) 

Using the LPA’ propagator (see eq. (11091) ) and the regulator of eq. (1201), making the 
change of variables p = p/n,v = q/n and cos 9 = p.q/p q, and performing the integrals over 
the remaining angular variables, one can write eq. © as: 


jf\^P)-~ 


K 


2vr— 


dvv 


d-1 


d -3 


X 


Z2(2vr) rf T(^±) (1 + ml) 2 J 0 Jo 

(2 — p + rjv 2 ) 


d6 sin 9(1 — cos 2 9) 2 


0(1 — v 2 — p 2 + 2vpcos9 ) + ( v 2 + p 2 — 2vpcos9)Q(v 2 + p 2 — 2vpcos9 — 1) + m 2 

(Bl) 


This expression is valid for arbitrary d, but we shall evaluate it only for d = 3. In order 
to take care of the 0 functions it is convenient to separate the calculation in two different 
regions: 2 < p and p < 2. In each case, one performs the 9 integral first, and then the 
integral over v. One gets: 
a) 2 < p 


jf’(K|p) = 


kZ|(2t)2(1+AJ) 2 


2 +f 


5 -2 

'3 +P 


1 

2 p 

1 
2 p 


1 + 4 + \ P + 


*1 + | + (p ~ 


'-bU’+s 





— (- 

k,Z'1(2tt) 2 (1 + ml) 2 \p 2 \3 


JL 

15 


+ 


p 4 


1 

15 


105 3 15 7 v ” 


(B2) 
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c) p < 2. 


4 3 \k-,p)-- 


„-i 


k ' J r] ijrh 1 2 K _ i 3 77 7rj fn 

: Z2(27r) 2 (l+m2)2 \ _1 + 4 + ^ + ^V2 _ 8 8~~ 


3?? -2 
-T p 


25// 

"48 


1 

+ 2 /j 
1 

2p 




1 - | - (p + 


v , v r~^2 

2 + 4 (P + V 


p + 1 + 


i ~\~(p- u-f+K^-^) 



K 


-1 


Z2(2vr)2(l+m2)2 ^(l + m 2 

ip 3 _ 


1-5 


7 P 2 


2 + 7/ - 2m 2 K + rjm 2 K 3 2(1 + 77 — 5?7i^ + r/m 2 ) 4 5 


6(1 + /Tl 2 ) 3 


5/ 3(1 + 77/2)2 

P 4 + 0 (p 5 6 7 8 9 10 11 ) 


15(1 +m2) 


(B3) 


ACKNOWLEDGMENTS 


Authors R. M-G and N. W are grateful for the hospitality of the ECT* in Trento where 
part of this work was carried out. 


[1] C.Wetterich, Phys. Lett., B301, 90 (1993). 

[2] U.Ellwanger, Z.Phys., C58, 619 (1993). 

[3] N. Tetradis and C. Wetterich, Nucl. Phys. B 422, 541 (1994). 

[4] T.R.Morris, Int. J. Mod. Phys., A9, 2411 (1994). 

[5] T.R.Morris, Phys. Lett. B329, 241 (1994). 

[6] J. Berges, N. Tetradis and C. Wetterich, Phys. Rept. 363, 223-386 (2002). 

[7] C. Bagnuls and C. Bervillier, Phys. Rept. 348, 91 (2001). 

[8] L. Canet and B. Delamotte, cond-matt/0412205 

[9] S. Weinberg, Phys. Rev. D8, 3497 (1973). 

[10] U. Ellwanger, Z. Phys., C62, 503 (1994); U. Ellwanger, M. Hirsch and A. Weber, Eur. Phys. 
J. C 1 (1998) 563; J. M. Pawlowski, D. F. Litim, S. Nedelko and L. von Srnekal, Phys. Rev. 
Lett. 93, 152002 (2004); J. Kato, arXiv:hep-th/0401068; C. S. Fischer and H. Gies, JHEP 
0410, 048 (2004). 

[11] U. Ellwanger and C. Wetterich, Nucl. Phys. B423, 137(1994). 


46 





























[12] For a review, see R. Alkofer and L. von Smekal, Phys. Rep. 353, 281 (2001). 

[13] T. R. Morris, Phys. Lett. B 334, 355 (1994). 

[14] G. Bayrn, J.-P. Blaizot, M. Holzmann, F. Laloe, and D. Vautherin, Phys. Rev. Lett. 83, 1703 
(1999). 

[15] G. Baym, J.-P. Blaizot, M. Holzmann, F. Laloe, and D. Vautherin, Eur. Phys. J. B24, 107 

( 2001 ). 

[16] J. O. Andersen, Rev. Mod. Phys. 76, 599 (2004). 

[17] G. Baym, J-P. Blaizot and J. Zinn-Justin, Europhys. Lett. 49, 150 (2000). 

[18] P. Arnold and G. Moore, Phys. Rev. Lett. 87, 120401 (2001). 

[19] V.A. Kashurnikov, N. V. Prokof’ev, and B. V. Svistunov, Phys. Rev. Lett. 87, 120402 (2001). 

[20] J. P. Blaizot, R. Mendez Galain and N. Wschebor, Non-Perturbative Renormalisation Group 
equations and momentum dependence of n-point functions (II), hep-th/0603163 

[21] J. P. Blaizot, R. Mendez Galain and N. Wschebor, Europhys. Lett., 72 (5), 705-711 (2005). 

[22] J. P. Blaizot, R. Mendez Galain and N. Wschebor, Phys. Lett. B632, 571-578(2006). 

[23] J. P. Blaizot, R. Mendez Galain and N. Wschebor, Non-Perturbative Renormalisation Group 

calculation of the self-energy of a scalar field, in preparation. 

[24] D.Litim, Phys. Lett. B486, 92 (2000); Phys. Rev. D64, 105007 (2001); Nucl. Phys. B631, 
128 (2002); Int.J.Mod.Phys. A16, 2081 (2001). 

[25] R.D.Ball, P.E.Haagensen, J.I.Latorre and E. Moreno, Phys. Lett., B347, 80 (1995). 

[26] J.Cornelias, Nucl. Phys., B509, 662 (1998). 

[27] L.Canet, B.Delamotte, D.Mouhanna and J.Vidal, Phys. Rev. D67, 065004 (2003). 

[28] J. Zinn-Justin, Quantum field theory and critical phenomena, Int. Ser. Monogr. Phys. 113, 1 

( 2002 ). 

[29] J. Polchinski, Nucl. Phys. B231, 269 (1984). 

[30] M. Bonini, M. D’Attanasio, G. Marchesini, Nucl.Phys. B409, 441 (1993). 

[31] L. Canet, B. Delamotte, D. Mouhanna and J. Vidal, Phys. Rev. B68, 064421 (2003). 

[32] T.R. Morris, M.D. Turner, Nucl.Phys. B509, 637-66 (1998); A. Ringwald, C. Wetterich, 
Nucl.Phys. B334, 506 (1990); C. Bervillier, hep-th/0501087 

[33] M. D’Attanasio and T. R. Morris, Phys. Lett. B409, 363-370 (1997). 

[34] N. Tetradis and D. F. Litim, Nucl. Phys. B464, 492-511 (1996). 

[35] M. Moshe and J. Zinn-Justin, Phys. Rept. 385, 69 (2003) 


47 




[36] R.Guida, J.Zinn-Justin, J. Phys. A31, 8103 (1998). 

[37] P. Arnold and B. Tomasik, Phys. Rev. A 62, 063604 (2000). 

[38] F. de Souza Cruz, M.B. Pinto, and R.O. Ramos, Phys. Rev. B64, 014515 (2001); Phys. Rev. 
A 65, 053613 (2002). 

[39] J.-L. Kneur, A. Neveu, M. B. Pinto, Phys. Rev. A69, 053624 (2004). 

[40] B. M. Kastening, Phys. Rev. A 69, 043613 (2004). 

[41] X. Sun, Phys. Rev. E67, 066702 (2003). 

[42] Note that depending on the choice of the regulator, not all fluctuations may be suppressed 
when k = A. However, for renormalisable theories, and if A is large enough, the effects of these 
remnant fluctuations can be absorbed into a redefinition of the parameters of the classical 
action. 

[43] The smoothness of n-point functions would be spoiled by a sharp cut-off regulator. The Litim 
regulator used in this paper would also lead to difficulties in high orders; however it allows for 
a Taylor expansion of n-point functions at leading orders in q 2 / k 2 , which is enough for our 
purposes. 


48 



